Heat conduction in one dimensional chains 
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We study numerically the thermal conductivity in several different one dimensional chains. We 
show that the phonon-lattice interaction is the main ingredient of the Fourier heat law. Our ar- 
gument provides a rather satisfactory explanation to all existing numerical results concerning this 
problem. 
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It is still an open and challenging problem to under- 
stand the macroscopic phenomena and their statistical 
properties in terms of the deterministic microscopic dy- 
namics. The crucial point is that how to connect the 
irreversibility with the time reversible deterministic mi- 
croscopic dynamics. One outstanding problem is that 
whether or not the heat conduction in one dimensional 
(1-D) chain obeys the Fourier heat law (normal thermal 
conductivity)? If yes, then under what condition. 

The first convincing result of the Fourier heat law in 
a classical system was given by Casati et al [1|. They 
studied the so-called ding-a-Hng model, which is a 1-D 
chain consists of the fixed equidistant hard-point parti- 
cle harmonic oscillators, and inbetween two fixed parti- 
cles there is a free particle. The particles have the same 
mass. The two ends of the chain are put into two thermal 
reservoirs. Classically, this system can be changed from 
integrable to fully chaotic by adjusting the system pa- 
rameter. They found that the key ingredient for the nor- 
mal thermal conductivity is chaos. Later on, Prosen and 
Robnik M have studied the ding-dong model by three 
different numerical methods and verified the Fourier heat 
law. The ding-dong model is a modification of the ding- 
a-ling model. The only difference is that in the ding-dong 
model the fixed harmonic oscillators are allowed to collide 
and there is no free hard-point particles inbetween. Fur- 
thermore, they have studied the temperature dependence 
of the thermal conductivity and found that it increases 
monotonically with the temperature. 

Most recently, Lepri et al M have studied the Fermi- 
Pasta-Ulam (FPU) (3 model. This model represents 
the simplest anharmonic approximation of a monoatomic 
solid. They put the Nose-Hoover thermostats act on the 
first and the last particle keeping constant temperature 
T-(_ and T_ , respectively. They shown that there exists a 
simple nontrivial scaling relation for the increasing num- 
ber of particles. The thermal conductivity, however, di- 
verges approximately as N^^'^. N is the number of parti- 
cles. They claimed that chaos is not sufficient to ensure 
the Fourier heat law. 

In this paper, we shall investigate the mechanism lead- 
ing to the Fourier heat law. In other words, we would 



like to answer the question: under what condition the 
heat conduction of a 1-D many-body Hamiltonian sys- 
tem having the Fourier heat law. To this end, we shall 
consider different models, such as the Frenkcl Kontorova 
(FK) model and the harmonic dissipative model etc.. We 
will show that by invoking a simple mechanism we can 
obtain a rather satisfactory explanation to all existing 
numerical results, qualitatively and quantitatively. The 
possible connection with the experimental results is also 
discussed. 

Normal thermal conductivity. — Either the ding-a- 
ling or the ding-dong model is more or less an artificial 
model. We would like to turn to a more realistic model, 
which is close to true physical system, i.e. the Frenkel 
Kontoroval model. It describes a particle (atom) chain 
connected by harmonic springs subject to an external 
sinusoidal potential. It has been widely used to model 
crystal dislocation, charged density wave, magnetic spi- 
rals and absorbed epitaxial monolayers etc. in condensed 
matter physics B . This model displays very rich interest- 
ing phenomena. However, we shall not discuss in detail 
all this properties in this paper, for more details please 
see Refs. HJq]. Our attentions are focused on the thermal 
conductivity in this paper. 

The existence of the thermal conductivity of this model 
has been proved by Gillan and HoUoway by using differ- 
ent numerical techniques 0]. The classical Hamiltonian 
of the standard FK model is 
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For convenience of numerical calculations, we shall scale 
this Hamiltonian into a dimcnsionlcss one, 
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By doing this, we have a new effective dimensionless tem- 
perature T. The real temperature T^ is related to T 
through the following relation. 
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where m is the mass of the particle, 7 the elastic con- 
stant, and b the period of external potential, which is 
unit after scaling; a the equilibrium distance of the par- 
ticle, it is /i = a/b after scaling; K = A/ja"^ is a rescaled 
strength of external potential. uJq — 7/m is the oscillator 
frequency, ks the Boltzman constant. In this paper, the 
winding number in FK model is kept at 1/3. 

It is helpful by establishing the above relationship (||). 
It can give us a very useful information about the cor- 
responding true temperature to that one we used and 
gain some physical insights. For instance, for the typical 
values of atoms. 
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we have Tr ~ 10^ — lO^T, which means that the room 
temperature corresponds to the dimensionless tempera- 
ture T about the order of 0.1 — 1. So, if T is very high 
such as to 10^, then the actural temperature is about 
10^ — 10^ degree, at this temperature the displacment of 
the particle from its equilibrium can be up to the order 
of 10, which we think is unrealistic for physical systems. 
Therefore, like Casati ei a/ [|l[, we always keep T at very 
small values in our numerical simulations. 

The Hamiltonian (g) is a very special case. In fact, we 
can write it into a general form, 

2 
H = J2h., H, = ?^ + Vix,^i,x,) + Uix,). (4) 
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Here, V{xi^i,Xi) stands for the interaction potential of 
the nearest-neighbor particles; U{xi) is a periodic exter- 
nal potential which is an analog of the lattice, and as we 
shall see later that it plays a crucial role in determining 
the behavior of the thermal conductivity. If U{xi) van- 
ishes and V{xi^i,Xi) takes the anharmornic form, Eq. 
(^ is then FPU /? model which has been discussed by 
Lepri et a^ [gj. Therefore, a variety of 1-D models can be 
put into the framwork of Eq. (Q) . By changing the form 
V{xi-i,Xi) and U{xi), we will obtain different thermal 
conductive behaviors. 

To study the heat conduction in 1-D model, we choose 
the same approach as that used by Lepri ei aZ |p[ , namely, 
two Nose-Hoover thermostats M are put on the first and 
last particle, keeping the temperature at r+ and T_, re- 
spectively. The equations of motion of these two particles 
are determined by. 



xn = —C-^N + In — Jn+i, 

c+ = 1^ - 1, c- = i^ - 1. 

The equation of motion for the central particles is. 



(5) 



where fi = —V'{xi^i — xi) — U'{xi) is the force acting 
on the particle, xq — Q and x^v+i = 0. 
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FIG. 1. Temperature profile for the FK model (H) with 
parameter K = ^. T+ = 0.3, T- — 0.2. The average is taken 
over 10® interval after the transient time 10*. The particle 
numbers are 300 (solid line), 200 (dashed line) and 100 (dot- 
ted line), respectively. 

We have carried out extensive numerical simulations 
with a large range of parameters of N, T± and K for 
a variety forms of T^(a;i_i,i) and U{xi). We used the 
seventh-order and eighth-order Runge-Kutta algorithm, 
which provides us much stable and more accurate re- 
sults than that from the usual fifth-order Runge-Kutta 
method. The spatial temperature profile for standard FK 
model is shown in Fig. 1 . It is clear that although our FK 
model has an additional external potential, if its strength 
K is sufficient large (compared with the temperature), we 
can obtain the same scaling relation as that obtained by 
Lepri et al. This scaling indicates that the temperature 
gradient scales as N~^. We have confirmed that this 
scaling relation is also true for many different modified 
FK models. For example, we have changed V{xi^i,Xi) 
to the anharmonic case as discussed by Lepri et al, or 
changed the external potential U{xi) to that one with 
higher harmonic term, such as. 
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The derivation of the heat flux of the ith particle dif- 
fers slightly from that of Lepri et al. The local heat flux 
J{x, t), which is defined by the continuity equation. Tak- 
ing the volume integration on both sides of this equation, 
we can obtain 



dV dV 
Ji - Ji-i = Xi — {xi,Xi+i) - ii-i-K {xi-i,Xi). (8) 

Thus the heat flux is defined by 
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Numerically, the time average J = {Ji{t)) is independent 
of the index i for long enough time. 
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FIG. 2. J X N versus the number of particles A'^ for dif- 
ferent models. T+ = 0.3 and T- = 0.2 for all cases. The solid 
circle represents the results of the FPU /3 model (/3 — 0.5). 
Open circle is the results of FK model given by Eq. (Q) with 
K = 5; solid square for FK model with an external potential 
(1^), with Ki = 5 and K2 — 15; and the open triangle for 
FK model with V = x^ /2 + (3x^/4:, U = -K cos{2iix)/{2-Kf , 
P — 0.9, K = 5. The lines are draw to guide the eyes. 

The N dependence of J x TV is plotted in Fig. 2 
for different models. As is easily seen that for differ- 
ent FK models (with different V(xi-i,Xi) and/or differ- 
ent U{xi)), as long as U{xi) is nonzero and at sufficient 
lower temperature, J x A'^ is a constant implying 1/J di- 
verges with N. Since the temperature gradient vanishes 
as N~^ as is shown in Fig. 1, thus the Fourier heat law 
is justified. 

Abnormal thermal conductivity. — Things become 
very different, if U{xi) vanishes. In this case, the heat 
conduction does not obey the Fourier heat law neither for 
the harmonic form V{xi-i^Xi) nor for the anharmonic 
form such as the FPU j3 model discussed by Lepri et al. 
Our results for FPU [3 model at very low temperature 
shown in Fig. 3 also demonstrate that J x N diverges as 
approximately N^''^ which means that the thermal con- 
ductivity diverges as N^^"^. This agrees with that of Lepri 
et al at much higher temperature. 

Based on the above results, we are convinced to con- 
clude that the key point of the normal thermal conductiv- 



ity is the periodic external potential, which is analogous 
to the lattice. 

If the lattice is absent, and the interparticle potential 
is harmonic, then no phonon-phonon interaction exists, 
thus the heat transfer would take place at the speed of 
sound and the thermal conductivity would be infinite as 
was pointed out by Debye in 1914. (However, if we add a 
dissipative term to the harmonic oscillator chain, then we 
could obtain the Fourier heat law, even though we have 
not lattice. This is because the dissipation, the heat ra- 
diates during the transport. Our numerical results have 
verified this. But we will discuss this more in detail in 
another paper j^].) 

In the case of having an anharmonic interparticle po- 
tential V{xi-i,Xi) such as that in the FPU model, the 
phonon-phonon interaction is produced due to the an- 
harmonicity. Although the temperature gradient can be 
formed, nevertheless, as is shown by the work of Lepri et 
al at high temperature as well as ours at low temperature, 
the thermal conductivity diverges. 

As long as the lattice exists, the phonons will be scat- 
tered by it and results in the thermal resistance, even- 
tually leads to the Fourier heat law. In the ding-a-ling 
model and the ding-dong models the fixed harmonic os- 
cillator plays the role of the lattice, whereas in the FK 
model, it is the periodic external potential. In these three 
cases the Fourier heat law is justified numerically. Thus 
we believe that it might be a general rule that if the 
phonon-lattice interaction is dominant, the heat conduc- 
tion will obey the Fourier heat law, no matter whether 
the interparticle interaction is harmonic or anharmonic. 

Temperature dependence of J . — As discussed above, 
the crucial point of the Fourier heat law is the phonon- 
lattice interaction. The mean free path of the phonons is 
determined by the density of lattice and does not change 
with the temperature. By increasing the temperature, 
more and more high energy phonons are excited, which 
results in the growth of the heat flux, thus the increment 
of the thermal conductivity. Whereas in the absence of 
lattice, increase temperature will produce more phonons, 
which in turn reduce the phonons' mean free path, con- 
sequently decreases the heat flux. Therefore, the tem- 
perature dependence behaviour for normal and abnor- 
mal thermal conductivity should be very different. Our 
numerical calculations exactly demonstrate this point. 

In Fig. 3, we plot the temperature dependence of heat 
ffux for different models. The particle number is kept 
at A^ = 100, and in all cases the temperature difference 
is fixed at AT = T+ — r_ = 0.1, thus J has the same 
behavior of the thermal conductivity k. For the FPU 
/3 model {(5 — 0.9, solid circle), the heat fiux decreases 
monotonically with temperature, whereas in the stan- 
dard FK model with {K — 5, solid triangle) increases 
with temperature. 



Another very important thing deserve noting is the 
case in which the anharmonicity and the external po- 
tential coexists. It seems that this case is closer to the 
real physical system than others. We have performed 
the numerical simulation by using V — x'^ /2 + I3x*/A and 
U = -Kcos{2TTx)/{2TTf in Eq.(|). The temperature de- 
pendence of J is shown in Fig. 3 (solid square) . The heat 
flow is affected not only by the phonon-lattice interac- 
tion, but also by the phonon-phonon interaction. At low 
temperature region, the factor determining the heat con- 
duction is the phonon-lattice interaction, therefore, the 
heat conduction obeys the Fourier heat law, but the heat 
flux is bigger than that case of the standard FK model 
(solid triangle) due to the anharmonicty which produces 
more phonons to transfer heat. The anharmonicity be- 
comes more and more important when the temperature 
is increased, this is why at higher temperature region, a 
relative flat region shows up in Fig. 3. Furthermore, it 
must be noted that our numerical results shown that for 
FK model shown in Fig. 3 the Fourier heat law is valid 
only at lower temperature region T < 1, at higher tem- 
perature the Fourier heat law broken down due to the 
reason mentioned before. 



ature region, whereas it decreases at high temperature. 
This can be understood well from the mechanism dis- 
cussed in this paper. 

In summary, by studying the dynamical equations of 
the 1-D particle chain, we understood more about the 
heat conduction mechanism. Our numerical results as 
well as others up to date confirm our conjecture that 
the phonon-lattice interaction is the key factor for the 
Fourier heat law. Only the phonon-phonon interaction 
cannot give rise to the Fourier heat law, instead we will 
have the abnormal thermal conductivity, i.e. the ther- 
mal conductivity diverges as the particle's number. In 
the former case, the thermal conductivity grows with the 
temperature monotonically, whereas in the latter case it 
decreases. 
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and the Hong Kong Baptist University Faculty Research 
Grant. 





• 
\ 


0.020 


-\ 




- \ 


0.015 


•x 




\ 


0.010 




0.005 


' ^/ 




.* 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

T 

FIG. 3. The temperature dependence of heat flux J for 
the FPU /3 model (/3 = 0.9, solid circle), the standard FK 
model Eq.(Q) (K=5, solid triangle), and the FK models with 
an anharmonic interparticle potential V(x) = x^ 12 -\- /3a;'*/4 
and external potential Kco's,(2'kx)I{2-k)'^ at (3 — 0.9, if — 5 
(solid square). The line is draw just for guiding the eyes. 

From many experimental results (see e.g. the book of 
Srivastava or other textbooks of solid state physics, 
such as that of Kittel ll^]), we observe that the thermal 
conductivity increases with temperature at lower temper- 
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